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We study photoinduced optical responses of one-dimensional strongly correlated electron 
systems. The optical conductivity spectra are calculated for the ground state and a photoexcited 
state in the one-dimensional Hubbard model at half filling by using the exact diagonalization 
method. It is found that, in the Mott insulator phase, the photoexcited state has large spectral 
weights including the Drude weight below the optical gap. As a consequence, the spectral weight 
above the optical gap is largely reduced. These results imply that a metallic state is induced 
by photoexcitation. Comparison between the photoexcited and hole-doped states shows that 
the photoexcitation is similar to chemical doping. 
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Since the discovery of high-T c superconductivity in the 
carrier-doped copper oxides, 1 doping effects in strongly 
correlated electron systems have been one of the most 
important subjects in the condensed matter physics. 2 In 
many cases, carrier doping is carried out chemically, for 
example, by the substitution of Sr for La in La2Cu04- 1 In 
addition to such chemical doping, photoirradiation is an- 
other way to introduce carriers into materials (photodop- 
ing). 3 Photodoping effects have been studied in many 
strongly correlated electron systems. 4-7 

One-dimensional (ID) halogen-bridged Ni complexes 
(Ni-X chains) are also known as strongly correlated elec- 
tron systems. 8 They are in the Mott insulator phase be- 
cause of large on-site Coulomb interaction on Ni sites. 
Recently, in Ni-X chains, photoinduced enhancement of 
a Drude-like low-energy component has been found just 
after the photoirradiation, suggesting the Mott transition 
by photodoping. 7 

It is well known that chemical doping of Mott insula- 
tors enhances the Drude weight in optical conductivity, 9 
and many theoretical studies on this effect have been per- 
formed, for example, on the basis of the numerical diag- 
onalization. 10-14 By contrast, the photodoping effect on 
optical responses in Mott insulators has not been studied 
so intensively. 15 Thus there remain several points un- 
solved. First, it is unclear whether photoexcited states 
are really metallic states accompanied with a completely 
different electronic structure from that of Mott insula- 
tors or not. It is established that Mott insulators have 
a characteristic electronic structure, the so-called upper 
Hubbard band and the lower Hubbard band. 16 Is such an 
electronic structure destroyed by photoexcitation? A sec- 
ond point is on the relation between photodoped systems 
and chemically doped systems. In both systems, a num- 
ber of carriers have significant influences on the physical 
properties. Thus we want to clarify how the photodoped 
systems are similar to chemically doped systems. 

To address these issues, we study the optical responses 



of the photoexcited state in the ID Hubbard model de- 
fined by 

l,a I 

where c, (q,<t) is the creation (annihilation) operator 

of an electron with spin a at site I, ni 1<7 — Ci Q^, and 
ni = m,i + ni t i. The parameter t denotes the nearest- 
neighbor transfer integral and U the on-site repulsion 
strength. In what follows, t is set to unity for simplicity. 
By using the exact diagonalization method, we calcu- 
late the optical conductivity and the Drude weight of 
the ground and photoexcited states. Our results show 
that the photoexcited state is metallic with a large Drude 
weight and a significant spectral weight below the optical 
gap A opt shifted from the holon-doublon continuum. As 
a consequence, the spectral weight above A opt is largely 
reduced. It is also found that the optical conductivity of 
the photodoped system is similar to that of a hole-doped 
system. 

The optical conductivity of the ground state j'i/'o) is 
given by 



where D is its Drude weight defined by 
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with K being the kinetic term of the Hamiltonian (1), j 
the current operator defined by j = i a [c { 
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cj ff c;+i )0 -), Eq the ground state energy, \ip n ) the n-th ex- 
cited state, and E„ the corresponding energy. The regu- 
lar component er rcs (cj) is defined by 
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where e gives a finite broadening and is set at 0.1 in our 
calculations below. 
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The optical conductivity of a photoexcited state is a 
key quantity. Here, we need to specify a photoexcited 
state. Experimentally, photoexcited states are created by 
the irradiation of a femtosecond-pulse laser, and they 
relax to the ground state within a few picoseconds. 7, 17 
Hence the true photoexcited states are non-equilibrium 
states, whose optical conductivity is generally difficult to 
be evaluated even numerically. Thus we here introduce 
some approximation on the photoexcited state: we regard 
the photoexcited state as the optically allowed, first ex- 
cited state IV'iopt}- In other words, IV'iopt) corresponds 
to the lowest peak of <7 reg (w). This assumption would be 
valid when we treat the system just after the irradiation 
of a pulse laser of energy equal to A opt . Then the optical 
conductivity of IV'iopt) with energy £i op t is given by 

(5) 



<7!(w) = Di<y(w) + a[ es (uj) + al es '(uj), 
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Here D\ is the Drude weight of IV'iopt) defined by 
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We treat the half-filled A-site chains and impose the peri- 
odic boundary condition for N~An and the anti-periodic 
boundary condition for N=An+2. In what follows, the 
system size is set at N = 12 unless otherwise noted. The 
Lanczos method is used for diagonalization. 

Figure 1 shows the optical conductivity <j(uj) and <j\ (uj) 
in the model (1). In <j(uj), the spectral weight distributed 
over a frequency range from approximately U — At to 
U + At is due to the single holon-doublon pair excita- 
tions 18 and called the holon-doublon continuum here- 
after. Below A opt , no spectral weight is observed except 
for the small negative D, which is visible only for [7 = 6. 
The appearance of a finite D is due to the finite-size 
effect, and D exponentially decreases with A. 19 

In (T\(lu), by contrast, several peaks appear below A opt . 
The most pronounced one is the Drude peak at oj = 0, 
implying the conducting property of the photoexcited 
state IV'iopt)- Another pronounced peak (we call it "A" 
below) is observed around u> ~ 1. The peak "A" is caused 
by the transition from IV'iopt) to the nearly degenerate 
state inherent in the ID Mott insulators. 20 Other small 
peaks are found below A opt [see the inset of Fig. 1(b) 
]. These are attributed to the transitions from IV'iopt); 
which is regarded as a holon-doublon pair excited state, 
to other holon-doublon pair excited states. The negative 
peak at uj — A opt corresponds to the transition from 
IV'iopt) to |Vo)- For w > A opt , the continuum is some- 
what similar to that in a(u>), but its weight is reduced. 
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Fig. 1. Optical conductivity in the ID Hubbard model, (a) for 
U = 6, and (b) for U = 12. The thin (thick) lines show a(u>) 
[ai(ui)]. The mid-gap region is magnified in the inset. 



In particular, the decrease of the peak height becomes 
significant as the peak energy is lowered, which is con- 
sistent with the experimental observation in the Ni-X 
chains. 7 

To clarify how the weight of the continuum ranging 
from U — At to U + At is reduced by photoexcitation, we 
calculate the integrated intensity defined by 



(uj)duj 



and 
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where we set uj\ and ui u such that these integrations only 
include the contribution from the holon-doublon contin- 
uum above A op t- We also calculate the total integrated 
intensity 
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and that of <7i(w), 



a(u))dio 



2N 



(Vol^lV-o) 



(11) 



roo 

li ot ^ ^ a 1 (u)du = — (V:,,,,,, 



I^IViopt). (12) 



These quantities are plotted in Fig. 2. We confirm 
that, in ui(uj), the integrated intensity from the holon- 
doublon continuum above A op t is largely reduced: the 
ratio I^ d /I hd is about 1/3 for U = 6. /} ld // hd increases 
with U, suggesting that the holon-doublon continuum 
becomes robust as U increases. It should be noted that 
jhdyjhd depends on the sys tem size [Fig. 2(c)]. 7 1 hd /7 hd 
increases with N owing to the decreasing density of 
photodoped carriers ppo Because IV'iopt) has two carri- 
ers (one holon and one doublon) , ppc is considered to be 
2/N. Hence the influence of photodoping becomes weak 
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as N increases. Note the large difference between /{ ld 
and l[ ot . In <7i(w), 7 1 hd is much smaller than l\ ot . Thus 
the low-energy components have dominant contribution 
to (7i(co). In a(ui), almost all the intensity is contributed 
from the holon-doublon continuum. 
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Fig. 2. (a) Integrated intensities (9)-(12), and (b) their ratios in 
the ID Hubbard model, as a function of U. (c) shows the size 
dependence. 



The large shift of the spectral weight from above A op t 
to low energies is the characteristic doping effect in the 
Mott insulator phase. 9 To confirm this in the photodop- 
ing case, we compare the results in the Mott insula- 
tor phase with those in the band insulator phase. For 
this purpose we hereafter employ the ID ionic Hubbard 
model (IHM) defined by 



Wihm = H + A^(-l)' 



ni, 



(13) 



where A(> 0) denotes a half of the level difference be- 
tween the neighboring orbitals. This model can describe 
both the Mott insulator phase (U > A) and the band in- 
sulator phase (U < A). 21,22 Figure 3 shows (t{uo) and 
<7i(w) in the Mott insulator phase (U = 12, A = 0) 
and in the band insulator phase (U = 0.1, A = 4). 
In G\{w) of the band insulator phase, the contribution 
from the hole-electron excitations 23 is similar to that in 



cr(u) except for the absence of the lowest peak around 
uj = 7.9 (= A opt ). This single difference between <j(lu) 
and ai(u>) comes from the addition of negative a[ eK '(u>) 
to a[° s (uj) in eq. (5). The bare contribution from the 
hole-electron pair excitation <7 1 es (o;) is almost the same 
as <t(oj), i.e., the peak positions and their heights are 
almost the same, suggesting that the band structure is 
almost unchanged by the photoexcitation. This is in con- 
trast to the situation in the Mott insulator phase. 
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Fig. 3. Optical conductivity in the ID IHM, (a) in the Mott in- 
sulator phase (U = 12, A = 0.0), and (b) in the band insulator 
phase (U = 0.1, A = 4). 



We have discussed the optical responses of the pho- 
toexcited state, and found that, in the Mott insulator 
phase, the photodoping substantially shifts a spectral 
weight from above A opt to low energies including the 
Drudc component. We now turn to a chemically doped 
system for comparison with the photodoped system and 
discuss their similarities and differences. The system 
doped with two holes is treated to be compared with 
the photodoped system with one holon and one doublon. 
We calculate the optical conductivity 



CT 2 hM EE D 2 hS{0j) + CT2h S ( W )» 



(14) 



where D 2 h is the Drude weight and cr^ 8 is the regular 
component of the ground state \1p2b) of the two- hole- 
doped system. The integrated intensities are also calcu- 
lated, which are defined by 
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and 



li°(Ih) = /. ^i(2h)(w)dw. (16) 



It is found that both doped systems have the reduced 
holon-doublon continuum above A opt [Fig. 4(a)]. The in- 
tegrated weight is comparable with FJ , and it ap- 
proaches I^ d as U increases [Fig. 4(b)]. In the low-energy 
part, there exists some difference. Although both sys- 
tems exhibit a large Drude component, the hole-doped 
system has no analogue of the peak "A" in ai(u>). In 
fact, D\ is smaller than D211 roughly by the amount of 
the peak "A" , so that the integrated low-energy weight 
I 1 ™ is comparable with li,° h w . 
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Fig. 4. (a) Optical conductivity of the undopcd (dashed line), 
photodoped (thin line), and 2-hole doped (thick line) ID Hub- 
bard model, (b) Integrated intensities (10), (15), and (16). 



To summarize, we have studied the photodoping ef- 
fect on the ID Mott insulator phase by using the exact 
diagonalization method. The optical conductivity shows 
that the photoexcited state in the Mott insulator phase 
has a large Drude weight, implying the metallic prop- 
erty. As a consequence, the contribution from the holon- 
doublon continuum above the optical gap, characterizing 
the ID Mott insulators, is largely reduced, suggesting 
the large deformation of the electronic structure in the 



Mott insulator phase. The comparison with the chemi- 
cally doped system shows the similarity between these 
doped systems. This fact would be modified if the op- 
positely charged carriers are attracted to each other by 
introduction of the nearest-neighbor repulsion. 
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